Typeset with jpsj2.cls <ver.l.2> 



Full Paper 



Quasi-one-dimensional spin-density-wave states in the presence of two 
commensurate potentials and an interchain misfit 

T. Itakura * 

Department of Physics, Kyoto-University, Kyoto 606-8502, Japan 
TEL.-075-753-3750 
(Received February 3, 2008) 

Spin density wave (SDW) states of a quasi-one-dimensional system with an incommensu- 
rate wave vector perpendicular to the chain have been studied in the presence of two kinds 
of commensurate potentials, which originate in a quarter-filled band and dimerization along 
the chain. In terms of a phase variable of the SDW order parameter, we treat classically 
the two-dimensional Hamiltonian, which includes both acoustic excitations with long wave 
length and a vortex excitation with short wave length. A phase diagram on the plane of 
temperature and chemical potential (where the latter corresponds to the deviation of the 
transverse wave vector from the commensurate one) exhibits a variety of states given by the 
commensurate SDW state without charge density, the commensurate SDW state with charge 
density, the incommensurate SDW state and the disordered state. 

KEYWORDS: spin density wave, phase, quarter-filled band, commensurate-incommensurate 
transition, renormalization group method, vortex, dimerization 

1. Introduction 

In several kinds of density wave systems, a commensurate- incommensurate transition has 
been studied theoretically by varying a misfit from the commensurate wave vector (or a 
chemical potential). Such a transition is followed by a formation of the soliton of the phase of 
density wave, since the energy gain of a commensurability competes with the increase of elastic 
energy driven by the misfit. Pokrovsky and Talapov have first studied such a phenomenon in 
a two-dimensional classical system by using a renormalization group method. 1 Such a method 
has been applied to the sine-Gordon model 2 ' 3 and the Coulomb gas model. 4 

In contrast to the one-dimensional case, it is crucial to take into account the effect of a 
vortex in the two-dimensional case, which appears at finite temperatures due to a thermal 
excitation. 5 The vortex destroys the quasi-long range order resulting in the disordered phase. 
The phase diagram of commensurate and incommensurate states in the presence of the vortex 
interaction has been studied by using the renormalization group method (for /x = 0) 6,7 and 
the exact result (for finite /i), 8 ' 9 where \x is the chemical potential corresponding to the misfit. 
The commensurability potential is characterized by the p-fold periodic potential with respect 
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to the phase. The case for p = 2 has been studied extensively by applying the renormalization 
group method to a one-dimensional quantum system, 10 in which the vortex and the misfit 
correspond to the anisotropic interaction of the spin freedom and an external magnetic field, 
respectively. 

These studies have been devoted to the system only with a single commensurate potential 
while that with two kinds of commensurate potentials is of interest to investigate the follow- 
ing quasi-one-dimensional conductors, Bechgaard salts. 11-13 The salts, which have a dimerized 
and quarter-filled conduction band, exhibit the spin density wave (SDW) state at low tem- 
peratures and undergoes a commensurate-incommensurate transition perpendicular to chain 
under effective pressures. 14 The interchain interaction including the misfit of the phase dif- 
ference between chains competes with intrachain commensurability energies where the misfit 
originates in the nesting of Fermi surface 15 and/or exchange energy due to the various kinds 
of interchain hopping. 4,16 When the repulsive interaction for the nearest-neighbor sites is 
large, the commensurability energy at quarter-filling leads to the coexistence of SDW with 
the charge ordering. 17 ' 18 However, the dimerization leads to an effectively half-filled band and 
is in favor of pure SDW state (i.e., the absence of charge ordering). In terms of a transfer 
integral method, the role of interchain misfit in the system with these two kinds of potentials 
has been examined for two-coupled chains 19 and for a quasi-one-dimensional system, where 
the phase diagram for the latter model has been obtained on the plane of an intrachain misfit 
and temperature. 20 

In the present study, we examine the effect of these two commensurate energies on the 
quasi-one-dimensional quarter-filled SDW system with a transverse misfit by adding dimer- 
ization and vortex to the previous model. 4 In §2, the classical two dimensional Hamiltonian 
for SDW state is derived by expanding the interchain coupling up to the second order with 
respect to phase difference between chains. The renormalization group equations are derived. 
In §3, we obtain the phase diagram of the commensurate SDW state and the incommensurate 
SDW state and the disordered state (induced by the vortex). In §4, discussion is given. 

2. Formulation 

We consider the SDW state of the quasi-one-dimensional electron system with the mean- 
field order parameter given by A^ cos(2/cf£o + Qi) where A^ and Oi denote the amplitude 
and the phase of the i-th chain, respectively. The effective Hamiltonian describing the charge 
fluctuation in the classical case is written as 




- W c cosA6i{xo) 



- W d wa.2Bi{x ) 



+ J cos(6»i(x ) - O i+ i(x ) + fi) 



(1) 
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where the ^4-term is an elastic energy while the W c - and W^-terms correspond to the com- 
mensurability energy at quarter-filling and that with dimerization, respectively. Equation (1) 
for two-chain has been studied previously. 19 In the case of = Jo = 0, the spatial vari- 
ation with a periodic array of (|,0, 1,0) is obtained for W c > 0, which is found for large 
nearest-neighbor repulsion. 17 ' 18 The Jo-term represents the interchain interaction with \x be- 
ing a misfit parameter corresponding to the phase difference between nearest-neighbor chains. 
We examine only the case of Jo < 0, because the replacement of 0i — ► 0% + tt in the even chain 
leads to the same result as that of Jo > 0. Since the calculation of eq. (1) is very complicated, 
we replace it by an effective one, which consists of acoustic excitation with long wave length 
and the excitation with short wave length. The former, which leads to the anisotropic two- 
dimensional SDW state at low temperatures, is obtained by expanding the Jo-term in eq. (1) 
with respect to 9i(xo) — 9i+i(xo) + \x up to the second order. The continuum limit is taken as 
&i(xo) d{xo,yo) with yo = ib where b is the lattice constant for interchain direction. The 
latter, which corresponds to the higher order terms and is dominant at high temperatures near 
the critical temperature, can be represented in terms of the vortex field, 6,7 which is relevant 
to the potential for 9 with a periodicity, 2tt. The potential for the vortex field is expressed as 

W v cos 2ttKc/)(x, y) , (2) 

which resembles a one-dimensional quantum system (see Appendix). In eq.(2), <p is the con- 
jugate phase variable defined as dd/dxo = i(3d</)/dyo and 89/dyo = —i/3~ 1 d(f>/dxo with 
(3 = ^(lJol/2^) 1 / 2 . The quantity K is given by K = y / 2A\J \/T , and the quantity W v is 
the coupling constant for the vortex field whose form depends on a microscopic model. The 
quantity Tq is the temperature. By use of dimensionless quantities, 




xq, y 



2/0 rp 

b ' 



J 



w c 



the effective Hamiltonian for eq.(l) is given by 4 
Wqid 



T 



- II 



dxdy 

90(x,y) 
dx 



+ 



J / d6(x, y) 



2T 



dy 



1 fw d 



— — cos 49(x, y) — — — — sin 29(x, y) 



T \W r 



T { W ' C ° S 27rK ^ x ' y ) 



1 (Wr. 



(3) 



For the convenience of applying a renormalization group method, eq. (3) is rewritten as 
Wqid 



n 



K 

2P 



dxdy 

dO{x,y) 
dx 



+ 



K(3 fd 6(x,y) 
dy 
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- 2j/ 4 cos 49(x, y) - 2y 2 sin 29(x, y) 



2y v cos 2nK(j){x } y) 



(4) 



where K = V2J/T, (3 = y/J/2, y 4 = 1/2T, y 2 = {W d /W c )/2T and y v = {W V /W C )/2T. The 
renormalization group equations are calculated by use of Coulomb gas method while the same 
results are obtained by the use of the method of Giamarchi and Schulz. 10 These equations are 
obtained as 

~ = -^M2p)-^J m+4**fi, (5) 
"' !!: ' - ^2 - -L") y 2 - 27ry 2 y 4 , (6) 

-2 



dl \ 7TK 

dy 4 

dl V" ttK 



V4~ny 2 , (7) 



^ = (2-nK)y v , (8) 
djl _ Airyl %ny\ 

Jl = ^ ~ -^ J i( 2 ^) - -^^i( 4 Ai) , (9) 

dln/3 47ry| 16vry| 

= -jf-Mty) + — ^^2(4m) (10) 

where yl = (3a 2 y 4 , y^ = (3a 2 y 2 , yZ = (3a 2 y v and 71 = /3a/x. The quantity a is short-length 
cutoff of the order of lattice constant where the constant is taken as unity. We note that the 
second order corrections to eqs. (6) and (7) are obtained by examining the response function 
up to third order, y 2 y4- 21,22 

3. Phase Diagram 

The renormalization group equations (eqs. (5)- (9)) are calculated numerically by choosing 
initial conditions as 1/K(= T/(2J) 1 / 2 ), 7i(= a^( J/2) 1 / 2 )), y 4 = a 2 K/A and y 2 = {W d /W c )y 4 . 
The vortex interaction term is taken as y v = ex.p(—7r 2 K/2) as found in the isotopic XY 
model. 6 In the present calculation, we take a = 0.5 and use an approximation for the Bessel 
function such that Jq(z) = 1, Ji(z) = z/2 for z< 1.8 and zero otherwise. These are possible 
choices to reproduce the well known result for y 2 = (see Fig. 1). 

From the solution of the renormalization equations, we obtain the following states. (i)The 
irrelevant 7i leads to two kinds of commensurate states. The relevant y 2 corresponds to the 
commensurate state of pure SDW (Ci-phase) while the irrelevant y^ and relevant yg corre- 
spond to the commensurate SDW state coexisting with charge order (Cn-phase) (ii) When 
71 is relevant, the incommensurate state (IC-phase) is obtained for the irrelevant yZ and the 
disordered state with the vortex (disorder-phase) is obtained for the relevant 77^. These fixed 
points are summarized in Table I where K* denotes a finite value. 

First, we examine the case that the dimerization is absent (y 2 = 0). In the present calcu- 
lation, we take a = 0.5 and use an approximation for the Bessel function such that Jq(z) = 1, 
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Table I. Fixed point values for respective phases. The parenthesis indicates the expected value. 

J\(z) = z/2 for z<1.8 and zero otherwise. These are possible choices, which reproduce well 
the previous results. 8,9 For small fx, the smaller a enhances the IC region while the larger a 
diminishes the IC region In Fig. 1, there are three fixed points corresponding to Cn, IC and 
disorder phases (Fig. 1). In the Cn-phase, the charge density becomes long range while, in the 
IC-phase, the gapless excitation due to a finite K*{= K(oo)) results in the power low decay of 
the response function of density waves. 4 In the disorder phase, the vortex, which is activated 
by thermal fluctuations, destroys the (quasi-)long range order and the response function for 
the density wave decays exponentially. There exists the incommensurate state between the 
commensurate state and the disordered state if fx ^ 0. The absence of the incommensurate 
state for ~p=0 is understood as follows. Equations (7) and (8) shows the irrelevant yl (y^) for 
1/K* > 7r/2 (1/K* < 7r/2) so that either yl or y^ is relevant for /7=0. We remark followings 
for other boundaries. Equation (8) shows that 1/K* is tt/2 at the boundary between disorder 
phase and IC-phase. The boundary between Cn and IC at T=0 (i.e., 1/K = 0) is given by 
(3li = 2V2/tt, which can be derived analytically by examining the condition for the forma- 
tion of the soliton. For > 2\[2/n, yl is always irrelevant, and the critical temperature 
for y^(l) — ► is almost independent of yl. The boundary at 1/K = is obtained from eqs. 
(6) and (8), which lead to d(yl 2 — 2lj2 2 )/dl = 4(yJ 2 — 2y2 2 ). When temperature is decreased 
(i.e., y 4 is increased) with fixed 2/2/2/4 = Wd/W c , it turns out that the Ci phase moves into 
the Cn phase for 2/2 < V2~V4 corresponding to the small dimerization. With increasing 1/K, 
the region of Q-phase increases, because the effect of 1/K on the scaling dimension for yl 
(i.e., 2 — 4/irK) is larger than that for 2/2 (he., 2 — 1/ttK). In addition, the high temperature 
corresponds to large 1/K due to 1/K = T/\J~2J so that Ci phase (Cn phase ) is obtained 
at high (low) temperatures. With increasing 1/K from the Cn-phase which is obtained for 
1/K* < tt/2, one always finds the Ci-phase, which is located in7r/2<l/K*< 2ir. 

Next we examine the state in the presence of both dimerization and quarter-filled com- 
mensurability energy, which exhibit four fixed points as shown in Table 1. Their states are 
given by two kinds of commensurate states. (Ci-phase and Cn-phase), incommensurate state 
(IC-phase) and disordered state (disorder-phase). In Fig.3, the phase diagram is shown on the 
plane of (3fi and 1/K where 2/2/2/4 = Wd/W c = 0.5. The Ci-phase (the commensurate state 
of the pure 2fcp SDW state) is found at high temperatures and for (3ji smaller than a criti- 
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cal value. The Cn-phase (the commensurate state coexisting with 4/cp CDW state) is found 
at low temperatures and the critical value of 0[i is larger than that of Ci-phase. Compared 
with Fig. 1, it turns out that dimerization has an effect of producing the region of Ci-phase 
around the triple point in Fig. 1, at which disorder, IC and Cn merge together. Note that the 
boundary between Cn and IC at T(= 1/K) = does not change within the visible scale even 
in the presence of both y<i and y v , which become irrelevant at low temperatures. With de- 
creasing temperature with fixed /3/U, one obtains following four kinds of successive transitions: 
(I) disorder -> Q -> C n , (II) disorder -> Q -> IC -»• Cn, (III) disorder -> IC -> C n , and 
(IV)disorder — > IC. However the tricritical point for Ci Cn and IC depends on the choice of a 
which leads to another transition of disorder — > IC — ► Ci — > Cn instead of (II). Therefore the 
analysis of such a detail is beyond the present scheme of the renormalization group equation. 
The remained three kinds of transition seem to be reasonable as the result of combining the 
commensurability of p = 2 and that of p = 4. 

4. Discussion 

We examined the effect of transverse misfit (chemical potential) on the quasi-one- 
dimensional SDW system which includes two-kinds of commensurate potentials coming from 
the dimerized quarter-filled band. The anisotropic two-dimensional SDW system is treated 
classically since the quantum fluctuations may have an effect of renormalizing the amplitude 
of W c and Wa but does not change qualitatively the form of the potential. We found the phase 
diagram, which exhibits a variety of states, on the plane of chemical potential and tempera- 
ture when the two kinds of commensurate potentials compete with each other. The transition 
from the disordered state to the incommensurate SDW state is the Kosterlitz-Thouless transi- 
tion. It is of interest whether such a transition in addition to commensurate-incommensurate 
transition is related to the sub-phases of (TMTSF) 2 PF 6 salt. 26 

Finally we note the phase diagram of Fig. 3, which is the result obtained in the classi- 
cal quasi-one-dimensional system with two-kinds of periodic potential competing each other, 
anisotropic misfit parameter, and vortex interaction induced by interchain coupling. This 
result in two-dimensional classical case has a common feature with the previous one, 20 (ex- 
cept for the state around the triple point) in which the interchain coupling is treated in the 
mean-field theory and the transfer integral method is applied to treat the thermal fluctuation 
exactly. In these cases, two kinds of potentials are the same but the misfit in the previous case 
is in the chain and the present case is perpendicular to the chain. Compared with the phase 
diagram (Fig. 8 in the previous case 20 ), we found the similarity in the sense that Ci (Ci) is 
the pure SDW state, Cn (C2) is the SDW state coexisting with the CDW and the disorder 
phase denotes the absence of the long range order, (i.e., (cos 9) = 0). The notable difference is 
between the incommensurate along the chain and that perpendicular to chain. It is considered 
that such a nice correspondence between these two results may be attributable to the proper 
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Appendix: Classical 2D Hamiltonian vs. quantum ID Hamiltonian 

Using the path integral methods, we show that the quantum one-dimensional (ID model 
is equivalent to the present classical two-dimensional (2D) model. Using the Stratonovich- 
Hubbard identity, 25 

2" 



/ 



V9i exp 



I 



VOiVTL 



x exp 



(A-l) 



the partition function for the two-dimensional Hamiltonian of eq. (1) is rewritten as 

z = f ve^vUi, 

x exp 



—— cos 46/j — — sin 26; 
T T 



+^cos(<9; - 6 i+ i + fj,) 
Jo 



d^ 

dx 



(A2) 



Next, we consider the one-dimensional quantum Hamiltonian with (Ax being the lattice con- 
stant) given by 



Tin 



zZ Ax 



In 2 

AA 3 



T 

gl/4 

2vr 2 a 2 



W d 

cos 49 j — — - sin 29 j 
T 



dx cos AB A 



fll/2 

2vr 2 a 2 



/ 



dx sin 29 



+ 



To 



+ jr cos(9j - 9 j+ i + fj) 



+ 



^ J dx Jocos(9i(x ) - 9 i+ i(xo) + 5) 



(A-3) 
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where [9j, Uj~] = iSjk and j denotes the lattice index along x-direction. Starting with eq. (A-3), 
the partition function for the Hamiltonian of eq. (A-3) is 25 

Z = Tr(e-P n «) 

i(To)\e- A ^H £ \0 j (n))(9 j (r 1 )\ 



= lim Tr(e 

N^oo 



- LIE 

-ArH q 



x e 



Mto)) 



- 1 



V9j VUj 



x exp 



To 



r r* 3 . 



Ax 



^Dj(T)-in,-(T) 



dr 



-^co S 4^(r)--^sin2%(r) 

^0 ^0 



+ 



Jo 



cos(^-(r) - Oj+i{r) + //) 



(A-4) 



where Ar = /3/iV and r is the imaginary time. Thus, the present classical Hamiltonian of 
eq. (1) is equivalent to the quantum Hamiltonian of eq. (A-3). By taking the continuum limit 
of H q along x-direction, namely 9j(r) — > 9 + (t,x) and — > f dx/Ax with x = jAx, one 
obtains following effective Hamiltonian, 



n = 



+ 



1 



— (d x e+ - v y + k p {d x t 



gl/4 

2vr 2 a 2 

Si; 



/■ 
/■ 



dx cos 
dx cos 9- 



Sl/2 



2vr 2 a 2 



/ 



dx sin 2# 



+ i 



2vr 2 a 2 . 

with the commutation relation, [8 + (x),d x 9-(x')} 

-///Ax, v p 



(A-5) 



-2iri5(x 



„/\23 



and if n 



T /(27T V /2^|J |), = —fi/Ax, v p = y/ | J 1/2 A, <? 1/4 = -2(7ra) 2 VF c /T and 5l/2 = 
2(7ra) 2 H / d/To. The Hamiltonian eq.(A-5) has been examined in the 1-D interacting quarter- 
filled electron system with dimerization. 22 

We comment on the interaction of the vortex in terms of 9—(x). In eq. (A-3), the Hamil- 
tonian is invariant under the transformation given by 9\ — ► 9i + 2n. The excitation expressing 
such a spatial variation corresponds to the short range topological excitation which is related 
to the vortex excitation. Actually, in continuum model, such an excitation is expressed as 
follows. Using the commutation relation, the vortex excitation is generated by unitary trans- 
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formation, in terms of V(x) = e ie ~( x \ that is V(x')9 + (x)V(x')^ = e ie -^9 + (x)e- ie - (x ^ = 
9+(x) + 7rsgn(x — x'). Then taking into account the symmetry of eq. (A-3), one can add the 
potential, J dxW v cos 9_{x) to eq. (A-5), since the term is invariant under the unitary trans- 
formation representing the lowest vortex excitation ( the vortex potential which has lowest 
scaling dimension). 24 Note that the correspondence between classical system and quantum 
system is given by x <-> r, y <-> x, b <-> Ax, 9 <-> 9+, 2irK(j) ^9^, K <-> 1/2ttK p and (3 <-> 
Thus we obtain the effective classical Hamiltonian eq. (3) which includes the vortex excitation. 
However, the microscopic calculation is needed to obtain the coefficient, W v . In the present 
paper, we put W v = ex.p(—Tr 2 K/2), which is the extrapolation deduced from the isotopic 
stiffness case. 6 
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